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ABSTRACT 


This report presents and interprets experimental data obtained from the 
Properties Measurement System (MPMS). Statics measurements yield the 
of-gravity of an unknown mass and dynamics measurements yield its 
matrix. Observations of the MPMS performance has lead us to speci ic 
criteria and an understanding of MPMS limitations. 


Mass 

center- 

inertia 

design 
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1. INTRODUCTION 

The Mass Properties Measurement System (MPMS), illustrated in Figure 1.1, 
consists of a serial kinematic mechanism with two intersecting revolute axes, z 0 
and zj, that intersect with fixed angle 04. The joint angles and rates denoted by 

0 1? 0 2 and 9j, 0 2 , respectively, turn about the respective joint axes. The angle 
rolls the table about z 0 and will be called tilt. 

The function of the MPMS is to measure the mass, center-of-gravity and second- 
order mass moments of an unknown mass placed on the turntable. The mass and 
center-of-gravity can be determined from statics measurements. The inertia 
tensor must be determined from the dynamics. Reference[l] provides a derivation 
of the MPMS dynamics equations. 

In all the experimental work reported here, we have fixed the angle 9 2 at some 

value during a measurement. In other words, 02 = 0 an< * = ® ^ or 
measurements performed in this study. 



Figure 1.1 The Mass Properties Measurement System 
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MPMS Dynamics 


The total torque x totall about the z 0 -axis equals 

x totali = x kei + x pei- ^ ^ 

According to ( 4 - 4 a) in [ 1 ], with = sin 04, x t = cos 04, 4 = cos 0 ^ s { = sin 0 i? 
i= 1,2, 14 = 0, r := r 2 and m := m 2 in that equation, the potential energy 

related torque equals 

t pe1 = 1110, g c {r, (c, c 2 - 1, s, s 2 ) + r y (-x, c 2 s, - c, s 2 ) + r* °l *1 >• (1-2) 
x equals the torque associated with the motion of the system (refer to ( 4 - 29 ) 

KEl 


in [1]), 


^kel 


£ T 


+ 1^01 s 2 (2 0 i 0 2 c 2 + s 2) 

• • • 

+ I yz (- 0 2 c 2 -2 0 1 c 2 x 1 + 0 2 2 s 2 + 20! 0 2 Xi s 2 ) 

+ (-®2 2 c 2 - 2 Q l c 2 X 1 - ®2 s 2 - 2 01 T 1 S 2) 

-1,^! (2 0 1 0 2 c 202 + 0i s 202 ) }• (1 " 3) 

The term I Tl describes the inertia of the first link and the other inertia terms, I xx , 
t t t t i describe the inertia of the second link plus the unknown 
mass. Ultimately, to determine the inertia tensor of the unknown mass will 
require knowledge of the inertia terms of the MPMS links. 

The analytical relation ( 1 - 2 ) provides the means for computing the center-of- 
gravity moment t g = mgr while ( 1 - 3 ) allows one to compute the inertia matrix 

I 2 foimed by the combined mass of link two and the unknown mass. 

MPMS Prototype 

The first axis of the MPMS prototype aligns with the diagonal of a unit cube 
(Fig Ure 1 . 3 ), hence, the angle the first axis makes to the floor equa s 


115 



Kennedy Space Center Aug 6, 1993 

MPMS... MEASUREMENTS NASA Faculty Fellow 


Keith L. Doty 
Kedron R. Wolcott 


atan2(V2,l)= 35.26°. The complement of this angle equals the twist oc^ = 
atan2(l,V2 ) ~ 54.74°. The body attached frame z-axis initially aligns with the 
axis of Oj. A rotation of the body by 0 2 = 120° will align the body attached frame 
x-axis to the first axis of the MPMS. Similarity, 8 2 = 240° will align the body 
attached frame y-axis to the first axis of the MPMS. This particular allignment 
permits determination of the diagonal entries of the inertia matrix by simply 
measuring the torque and position of the first axis. Computation of the off 
diagonal terms, however, requires 0 2 to be time varying and measured. 



Figure 1.3 MPMS prototype with aj = atan2(V2, 1) « 54.74°. 


2. COMPUTATION OF THE CENTER-OF-GRAVITY 

x 

To compute the center-of-gravity vector r = m g , the total mass m and 

gravitational constant g must be known. In the experiments to follow, we treat 
Oj m g r x , cfj ra g r y and Oj m g r z as the unknown terms and fix both axes. 

Measuring the reaction torque on the first axis for different configurations we 
can determine the actual values of r x , r y and r z when the solutions are divided by 

Oj m g. Details of this procedure follow. 
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MPMS Statics Coefficient Matrix 

The coefficients ofo^g r x , Gf m g r y and a A m g r z m (1-2) equal 

c x := Ci c 2 -t iSi s 2 , c y := -X! c 2 Sl - c t s 2 and c z :=r z a lSl , (2-1) 

respectively. At 0 2 = jjl, \i + 120° and p + 240°, where p equals an arbitrary, but 
fixed, angle, the coefficients form the coefficient matrix 


Cf = 


r 

- c i ( c p + s p 

V C 1 (" C P 


c p c l' x i s p S J_ 

)/2- (V3 - s^) sj/2 

+ y[3 s^/2 + Xi (V3 c^ + s^) s ! /2 


-C 1 s p ‘ c p x i S 1 


-ci (V3 c^ - s^/2 + 


x, (c^ + V3 s u ) sj/2 


<*1 S! 
Gi S! 
G i SI 


A 


J 


( 2 - 2 ) 


Cl (V3 C^ + s^)^ + x t (cjjl - V 3 Sp.) si/2 

A note about these choices for 0 2 . Physically, the angles 0 2 = 0, 120° and 240° 
were chosen to align axis one with three orthogonal axis F b selected in e 
unknown mass. These choices of 0 2 permit calculating the diagonal entries of the 
inertia matrix, expressed in F„ coordinates, from torque and positio 
measurements on the first axis alone. Serendipitously, these “gles allows us to 
find an MPMS configuration for which the condiuon number of C f eq 
fact any three angles separated by 120° will allow the MPMS to be configured 
such that a condition number of one can be achieved. Recall that a mamx C f w.th 
a condition number of one guarantees the optimum numerical stability when 
inverting Cf. 

By measuring the torque at joint one for each of the three angles selected for 
0 2 , we obtain the linear system of equations in die center-of-gravity u ow 


(2-3) 


(TlO 



T 21 1 

= m g Cf 

r y 

v T 3iy 


Vz) 


Solving (2-3) produces 
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ly = 


c l m g 


Cf 1 T 21 


(2-4) 


Two central issues arise: 1) For what configurations is C f invertible and 2) how 
well conditioned is C f for those situations where it is invertible? 

The coefficient matrix C f is invertible as long as 

3/2 

det[Cf] = TT (6 + 2 c 2ttl - c 2(ar 9l ) + 2 c^,- c 2( a 1+ 9 ,)) H * °- < 2 ’ 5a > 

Notice that the determinant does not depend upon the initial angle p. selected for 
e 2 as long as the angles differ by 120°. The determinant of the coefficient matrix 

can also be written simply as 


det[C f ]=— (l-a/s 1 z )a 1 s 1 


(2-5b) 


by using the not so obvious, but directly verifiable identity 

6 + 2 c 2(Xl - c 2(otr 01 ) + 2 c 29i - c 2(ttl+ 6l ) = 8(1-0! sj ). (2-6) 

The condition o lSl *0 means that the angle 0, and the twist a, between the two 
axes cannot be zero or a multiple of it. These design constraints make physical 
sense. Positioning the turntable at three different angles with Oi - 0 Wl11 not 
change the z-component of the center-of-mass, and, therefore, cannot provide 
information about r r The same observation applies to the case where the twis 

a.| equals zero. 

The term (1 - O , 2 S ! 2 ) = 0 when «,= 0,= ± \ it. The MPMS with this physical 
configuration can only measure r z . 


Design Criterion 

The MPMS configuration and structure must satisfy Sj ^ 1 and Sj ^ 0. 
In particular, the MPMS must have a non-zero twist a j between the two axes. 
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In the course of experimentation, sometimes the computed solutions based on the 
data vS wUdly Analysis of the causes revealed that one error source was the 
ill-conditioning of the numerical inversion of C f . 

Condition Number of the MPMS Statics Coefficient Matrix 

The condition number c n for C f equals 




^max 


l min | 


(2-7) 


where and X min are the maximum and minimum eigenvalues, respectively, 

of C x C f We now calculate the condition number c„ of C f as a function of the 
twista, between the shafts and 0„ the angle of the first shaft. This will allow us 
to determine the angle settings for which 

Using (2-2), direct computation yields 

( Aq 


Cj C f = 


V 


0 

o x 2 
0 o 


0 'N 
0 

^3 ) 


where X { , X 2 , X 3 clearly equal eigenvalues of C] C f and compute to 

X { = ( 6 + 2 c 2ai - c 2(0C r 0!) + 2 C 20! - C 2(«i+ 0l> 

X 2 = ^ (6 + 2 c 20tl - c 2(0Cl _ 01 ) + 2 c 201 - c 2(otl+ 01 ), 


X 3 = 3 of Sj 2 , 


(2-8a) 


(2-8b) 

(2-8c) 


and do not depend upon the arbitrary angle p added to 0°, 120° and 240° to 0 2 . 


Define y = 


V 


X- 


Applying (2-6) produces the simple expression 
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y:= 



1 -<M 2 »i 2 


(2-9) 


Whenever X 3 < Xj, the condition number c n of C { equals y. On the other hand, 
if Xj < Xj , the condition number c n of Cf equals l/y, 

if -y < 1, then c n = 1 ly else c n = y. (2-10) 

Selecting the Twist and Tilt to Optimize Condition Number 

A plot of c n versus 0! appears in Figure 2.1 for a, = atan2(l,V2 ) - 54.74°, the 
angle employed in the prototype (Figure 1.3). Details of the W-shaped port, on 
of the curve is illustrated in Figure 2.2. 


= atan2(l,3/2) 
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04 = atan 2 (V 2 , 1 ) 



(radians) 

Figure 2.2 Details of the condition number 04 = atan 2 (V 2 , 1 ). 

The condition number c n = 1, its optimum value, at that value of for which 
c n - 1 = 0. This occurs at 0 lopt = 0.7853981633974484 180/n = 45°. 

A plot of c n versus 0 t appears in Figure 2.3 for 04 = tc/ 2. Details of one period 
of the curve is illustrated in Figure 2.4. 

1 


C n 04 = 2 * 



0 jl (radians) 

Figure 2.3 The condition number c n of the coefficient matrix C f vs the angle 0! 
of the first shaft for 04 = 5 n - 
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1 

c n <*1 = 2 n 



The condition number c n =1 at that value of for which c n - 1 = 0. This occurs 
at 

0 lopt = 0.6154797086703873 180/71 = 35.2644 = = atan2(V2, 1 ). 

MPMS Configurations with Condition Number of One 

A general plot of the condition number (Figure 2.5) illustrates the nature of the 
condition number when the three angles selected for 0 2 split the circle into three 
120° arcs. The plot has been truncated at c n = 5 to display the important 
characteristics of the function. 
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Figure 2.5 Plot of the condition number, clipped at c n -5, as a function of 0; and 
the angle a, between the first and second axis. 

The condition number plotted in Figure 2.5 applies for any three angles chosen 
for 09 which are separated by 120“ The funct.on c n is symmetric in 0, and a, 
and has period of it along each axis. Figure 2.6 illustrates a single quadrant of 
Figure 2.5 in more detail and Figure 2.7 depicts an overhead view. 
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Condition Nu/nber 



Figure 2.6 Plot of one quadrant of the condition number clipped at c n 5. 

When the condition number is one, all the eigenvalues of C f must equal one. 
In particular, for c n = 1, (2-9) implies 

7t 3 = 3o 1 2 s 1 2 =3 sin 2 (oti) sin 2 ^) = 1. (2-11) 

Furthermore, if X 3 = 1, (2-9) and (2-8) easily proves that X l =X 2 = \ and, 
hence, c n = 1. Consequently, (2-11) provides a design equation for the MPMS . 
For example, the minimum tilt possible for c n = 1 occurs when a x 2 = 1, »- e -> 
(Xj = ±90°. The minimum angle 0j must, therefore satisfy s x 2 = 3 , yielding = 
±35.26°. 



Kennedy Space Center 
MPMS...MEASUREMENTS 


Aug 6, 1993 
NASA Faculty Fellow 


Keith L. Doty 
Kedron R. Wolcott 


15 


The crease at the bottom of the well in Figure 2.7 
which the condition number of the coefficient matrix 


depicts the curve (2-11) on 
equals one. 






nn 


Condition Humber 


Alpha 1 



Theta 1 


Figure 2 6 Top view of one quadrant of the condition number clipped at 
c„=5. 

3. STATICS MEASUREMENTS 

Another source of error discovered in M P MS cem of.grav,,y measurements 
^ armracv of the torque sensor attached to zq. 4 

was the absolu J approximately 1900 in-lb. Measured torques 

MSb 1 " 1“ r S ^ -f ^ — number 

° f ^^“'^0,^2?. the error in the values^ the center-of- 
gravity components could be matrix^ 

Sr.wtt S — 1 °a°n g les that would yield a condition number of one. 
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Statics Measurement Criterion 

The three turntable angles (9 2 > should trisect the circle and the tilt angle 9, 
should be selected to satisfy (2-11) for a given a 1 . 


Tilt Axis Torsion 

The dynamics measurements later will require oscillating the tilt axis at some 
prescribed frequency. Any significant torsion spring action of the tilt axis would 
affect the dynamics measurement. Essentially, we wanted to determine the M 
natural frequency to ascertain whether the tilt axis torsion would have to be 
modeled in the dynamics. The prototype MPMS, upon which all experimen s 
were performed, possesses the twist angle ccj = 54.73°. 



Figure 3.1 depicts the relationship between the angular deflection 60 i of the tilt 
axis as a function of the torque Tj exerted about that axis. The experimental 
compliance c = 1/k of the tilt axis equals c = 0.0019°/in-lb. Assuming frictionless 
operation, the natural frequency of the MPMS with the second joint fixed 
approximately equals 16 Hz for a 50 lb force, 4 in from the first axis, 
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32.2*12 

0.0019*50*43 


= 16 Hz. 


(3-1) 


In order to ignore the tilt axis torsion spring effect in the dynamics 
measurements, we should oscillate the table well below the resonant fr ^ uen ^ ° 
the second link and payload. A payload 9 times larger than the 50 lb fo ce, 
however would produce a resonant frequency of about 5.3 Hz. Such a lo 
resonance would probably require modelling the shaft torsion for the dynamics 

measurements. 

Center-of-Gravity Measurements 

The minimum-norm, least-squares values for the 

turntable plus an unknown mass placed on the turntable, yielded a first-mass 
moment v ec tor of m g r = I 160 -41 284 ] (in-lb). The vector m g r was 

calculated from the average values of 3000 torque samples per ex P« r, ” ent ' 
from 30 experiments. The 30 experiments consisted of permutations of 10 
measurements taken at each of the three MPMS configurations {( 0, , 0 2 ) I - I 
(35°, 0) , (-35°, 120°), (-35°, 240°) ). 

Since the torque sensor used in the prototype MPMS generated absolute errors of 
19 in-lbs, small, first-order mass moments measurement were exceptionally error 
prone. Figure 3.2 illustrates this point quite clearly. The three curves indicme . he 
error of each experiment about the solution m g r = [ 160 -41 284] (l )■ 

The "blue” plot (the one that starts in die upper left-hand corner of the diagram) 
indicates the most error in measurement. Although 19 in-lb almost equal 50% of 
the measured y-axis first mass-moment of -41 in-lb, the huge number of samples 
over which averages were computed produced consistent resu wi ° 

lamer values of the x-axis and z-axis first mass-moments produced more accurate 
roSgs within ±1% for m g r x (the "red” plot has a negative relative error for 

the first experiment) and within ±0.5% for m g r z ("green’ plot). These results 
are consistent with the fact that the absolute error of 19 in-lb has much less effect 
on the larger torque measurements. 
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Figure 3.2 


Experiment Number 

Percent error about mean of first order mass moment. At the far left 
the top plot is "blue", the middle plot "green" and the lower plot red 


4. DYNAMICS MEASUREMENTS 

The dynamics equation (1-3) for 04 = 54.73° (cq = \ 3 » T 1 = V 3 )’®2 _ “ 

0 , produces 

X kel = {Iri+Im} 0l. (4 " 1} 

where 


I m ■=\(2l xx s 2 2 + 2 ly f C 2 2 + l zz ~ 2 I xy S 292 - 2 3/2 I xz S 2 - 2 3/2 Iy Z c 2 ) } 0i .(4-2) 


To test the MPMS prototype, we measured the angular position and total 
torque x t0tall = x KE1 + x PE1 at joint one versus time in order to calculate the total 

inertia seen by the first joint, 
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{I T i + I m } 


TWall X PE1 


(4-3) 


e 


i _ t, v o factor from 10 to 60 , hence, the torque 

SET, TXTmt two nearly e q na. numbers 

numerical errors. ^ ^ 

significant problem^ To ^ second link and the 

the gravity torque x PE1 ~ 0. We > then « nanu y = ^ measured the reS ultant 

unknown mass pomt at about 8.7 Hz about the PE1 

torque and angular position. 

We compute the angular acceleration 0, from 6, by the central difference 

formula [1], + _ 30 0 , (t 0 ) + 16 9 ,(t.i) - + o(h 4 ).(4-4) 

.0,(t o ) =— 12 h 1 

where h = t, - t M = 4 ms and 0(h 4 ) indicates errors on the order of h 4 . 

Figures 4.1 and 4.2^depict semm^ed irTthe 

promtype^esulted Indarge Ltantaneous errors when measuring the small inert.. 
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Position vs. Time 



Figure 4.1 A plot of 0j vs time. 



Figure 4.2 A plot of x totall vs time. 
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torques. The mean of the absolute value of the torque lx kel l and position derived 

angular acceleration 1*0 x | were computed for each experiment consisting of 
5000 points. The ratio of these values were then taken as the estimate or t e 
measured inertia {I T i + I m 1 » 


{ I T i + f m ) 


It, 


kel 


(4-5) 


fe 


Figure 4.3 indicates the inertia measurement percent error for 10 different 
experiments with respect to their mean value. The error is bounded ±10% about 
a mean of 4.34 in-lb-s 2 These large error bounds can plausibly be attributed to 
the torque sensor sensitivity. The measured torque values in 
than 10% full scale, hence, the torque sensor, which is rated 1% accu y 
scale, produces approximately 10% error in the measurements. 
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5. OBSERVATIONS OF MPMS PERFORMANCE 

With single torque sensor and position sensor, one can measure die center of 
gravity and the moments of inertia about three orthogonal axes. However, one 
must minimally move and measure the velocity and acceleration of the second 
axes in order to obtain the full inertia tensor. These facts were established in [ ]. 

To obtain accurate center-of-gravity measurements one should select MPMS 
configurations where the condition number c n of the coefficient matrix equals 
one In other words, the choice of angles selected for measuring the center-of- 
gravity matter significantly. A poor choice of angles can multiply 
error by a factor 1 < c n < ~. This explained why, in some of our earli 

experiments the errors where 600%. 

Torque readings for both statics and dynamics measurements at the lower hrrut of 
the torque sensor resolution also had significant impact on the measurem 
reliability. The location of the center-of-gravity of the unknown mass affects th 
size of the torque readings. This suggests placing the mass center away from the 
z t -axis in order to drive the force sensor harder and produce more accura 
results Taking these details into consideration center-of-gravity measurements 
we made highly repeatable measurements (±0.4%) for a specific configur 

and setup. 

Gravity torque dominates dynamics measurements and extreme measures must be 
taken^o mLnize its affects. Dynamics measurements were made at a balance 
point of the turntable using “high frequency oscillations ( ' 

combination of the smaller gravity terms about th ' baan “ 
increased inertia effects for high speed oscillations enabled us to more accmate^ 
measure the much smaller inertial terms than a P nce ,^ ar 

Higher frequency oscillation are precluded because o Wor a eive” 
16 Hz Dynamics measurements were repeatable (_3. ) § 

configuration and setup and stayed within ±10.0%, even though die measurements 
were near the lower resolution limit of the torque sensor. 

Stiction and friction also appear to play an important role in inertia 
measurements on tire MPMS. Stiction was particularly 

the turntable. These effects were not modeled in this effort, but should be tak 
into account in a later prototype. 
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6. RECOMMENDATIONS 

In [1] it was recommended to make the twist 04 = 90°. Further support for this 
design change arises in the condition number c n analysis provided earlier. That 
analysis proved that the minimum tilt angle of 35.26° for c n = 1 can on y e 
achieved when 04 = 90°. By comparison 04 = 5.74° in the current prototype an 

the tilt angle must be 45° in order for c n = 1. 


At this point the limitations of the current prototype appear to be well understood 
and further experimentation is contraindicated. If management decides 
continue with the MPMS concept, the follow-on prototype will require some 
modificatins to correct the discovered problems with the first model. 


The prototype MPMS permitted only manual actuation of joint one. A second 
generation MPMS should : 


1 . Include automatic operation of both axes, 

2. Match torque sensor with expected loads, 

3. Utilize low friction bearings, 

4. Add a position sensor on the second axis, 

5. Incorporate a torque sensor on the second joint, 

6. Explore concept with potential customers 


7. CONCLUSIONS 

Theoretical results verify the MPMS concept in principle, but experimental 
results indicate several engineering limitations in the current ^P^mentatiom n 
particular the torque sensor must be able to accurately measure a broad range o 
values. The gravity terms dominate the dynamics and this fact makes it difficu 
to obtain good inertial measurements. 

To obtain the optimum performance from the MPMS when measuring center-of- 
gravity: 

1. Place the unknown mass's center-of-gravity away from the turntable 
center. 


133 



Kennedy Space Center Aug 6, 1993 Keith L. Doty 

MFMS... MEASUREMENTS NASA Faculty Fellow Kedron R. Wolcott 


2. Measure the center-of-gravity every 120° turn of the turntable. 

3. Tilt the turntable to the optimum angle to produce a condition number of 
one for the coefficient matrix. 

When measuring inertia : 

4. Place the unknown mass so that the turntable axis passes through its 
center-of-gravity. 

5. Oscillate the turntable well below the resonant frequency of the MPMS 
about the null gravity torque configuration . 

Items 1 and 4 clearly compete with one another. Placement issues have not been 
addressed here and would need further exploration should the MPMS project 
proceed. 
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